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We construct the external metric of a slowly rotating, tidally deformed material body in general 
relativity. The tidal forces acting on the body are assumed to be weak and to vary slowly with time, 
and the metric is obtained as a perturbation of a background metric that describes the external 
geometry of an isolated, slowly rotating body. The tidal environment is generic and characterized 
by two symmetric-tracefree tidal moments £ a b and B a b, and the body is characterized by its mass M, 
its radius R, and a dimensionless angular-momentum vector 1. The perturbation accounts for 

all couplings between \ a and the tidal moments. The body’s gravitational response to the applied 
tidal field is measured in part by the familiar gravitational Love numbers and A™ ag , but we find 
that the coupling between the body’s rotation and the tidal environment requires the introduction 
of four new quantities, which we designate as rotational-tidal Love numbers. All these Love numbers 
are gauge invariant in the usual sense of perturbation theory, and all vanish when the body is a 
black hole. 


I. INTRODUCTION AND SUMMARY 

In an earlier work (Q> hereafter referred to as Paper 0), the geometry of a slowly rotating black hole deformed 
by tidal forces was determined under the assumption that the tidal forces are weak and vary slowly with time. 
The construction relied on perturbative techniques, and the metric of the deformed black hole was expressed as a 
perturbation of a background Kerr metric with mass M and dimensionless spin <C 1. The perturbation introduces 
new terms to the metric, and these are constructed from tidal moments £ a b and B a b that provide, at leading order 
in the tidal deformation, a complete characterization of a generic tidal environment. In particular, the perturbation 
includes terms that arise from the coupling between \ a and £ a b, and the coupling between y a and B a b', these capture 
all consequences of the dragging of inertial frames on the black hole’s tidal deformation. The metric was cast in light- 
cone coordinates (v,r,Q,<j>) that possess a clear geometrical meaning: the advanced-time coordinate v is constant on 
light cones that converge toward the black hole, the angular coordinates ( 6 , <j>) are constant on the null generators of 
each light cone, and the radial coordinate r is an affine parameter on each generator. 

In this paper (paper I) we replace the black hole by a compact body of mass M , radius R , and dimensionless spin 
C 1, and rely on the fact that to first order in y“, the external geometries of the unperturbed bodies are identical 
(differences occur at second and higher orders). We generalize the perturbative solution of Paper 0 to account for the 
presence of matter inside r = R. The new solution is no longer required to be regular at r = 2 M (which marks the 
position of the black hole’s horizon), and as a consequence we find that it contains terms that were absent from the 
black-hole metric. These describe the body’s response to the tidal deformation, and they come with dimensionless 
multiplicative constants known as gravitational Love numbers. In the case of a nonrotating body (y“ = 0) M , two 
types of Love numbers appear in the metricQ a gravito-electric Love number K|* that measures the body’s response 
to the tidal forces associated with £ a b , and a gravito-magnetic Love number I\ ™ ag that measures the body’s response 
to those associated with B a b . With rotation included to first order, we find here that the description of the external 
geometry of a tidally deformed body requires the introduction of four new numbers, which we call rotational-tidal 
Love numbers. The first two, denoted £ q and 3° below, are associated with terms that couple \ a to £ a b in the metric; 
the remaining two, denoted *8 q and j?°, are associated with terms that couple \ a to B a b ■ The label q indicates that 
(£ q and © q are associated with quadrupolar terms in the metric; the label o indicates that and M° are associated 
with octupolar terms. All six Love numbers are gauge-invariant in the usual sense of perturbation theory: While 
the form of the solution may change under a gauge transformation (understood as an infinitesimal transformation 
of the background coordinates), the value of the Love numbers are preserved. And our construction guarantees that 
A'f 1 = A'“ ag = 0: q = = 55 q = M° = 0 for a slowly rotating black hole. 


1 A different set of Love numbers, klf and fc™ ag , was introduced in Ref. [)» . It is related to the set used here by kkf = (2 M/ Rf' and 
fc™ ag = (2M / Ti f' R)" as . The former notation for the gravito-magnetic Love number was unfortunate, because Favata |5| has shown that 
Klf ag scales as ( R/M ) 4 instead of (R/M) 5 . 
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Our considerations in this paper are limited to the exterior geometry of a slowly rotating, tidally deformed body. In 
this context, the new Love numbers (2: q , S' 0 , 93 q , .ft 0 ), like the old Love numbers (Iff 1 , IL™ as ), must be left undetermined. 
The determination of the Love numbers requires the construction of an internal geometry based on an assumed model 
for the body’s matter content, and a match of the internal and external geometries at the matter’s boundary; this 
procedure will be detailed in Paper 10. The numerical values adopted by the Love numbers reflect the details of 
internal structure, and an external measurement of the tidal properties of the body can therefore reveal otherwise 
inaccessible aspects of the interior. 

This observation has motivated a recent surge of activity in the development of a relativistic theory of tidal deforma¬ 
tion and dynamics, in the context of the measurement of tidal effects in gravitational waves emitted by neutron-star 
binaries |a-ll9| and during the capture of solar-mass compact bodies by supermassive black holes |20h25| . The Love 
number iff 1 of a neutron star was implicated, along with its moment of inertial I and its rotational quadrupole 
moment Q , in the remarkable J-Love-Q relations [26h32| . and tidal invariants have been inserted within point-particle 
actions to account for the tidal response of an extended body [tMitI . The geometry of a tidally deformed black hole 
was constructed beyond the small rotation limit |3S, and beyond perturbation theory j4(j. 

We begin our work in Sec. [IT]with the construction of tidal potentials obtained from Xa , £ab, and B a b- These form 
the fundamental building blocks of the perturbed metric, which is calculated in Sec. IIIII in the light-cone coordinates 
(v,r,d,(j)). To illustrate the features of this metric, it is helpful to introduce Cartesian coordinates x a = ( x,y,z ) 
defined in the usual way from the spherical polar coordinates (r, 9 , <j >). In this notation, the time-time component of 
the metric obtained in Sec. mu takes the form of 
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in which e a bc is the completely antisymmetric permutation symbol, all indices are raised and lowered with the Euclidean 
metric 5 a b, the notation ( abc) indicates symmetrization of all indices and removal of all traces, and ellipses designate 
relativistic corrections of order 2 M/r and higher. 

The first two terms in g vv originate from the background metric, and the following terms proportional to £ a b x a x b 
describe a quadrupolar tidal deformation of the geometry; the decaying piece involving ATf 1 represents the body’s 
response to the applied tidal field. The next sequence of terms describes another quadrupolar deformation that arises 
from the coupling between x a and £ a b ; the decaying piece involving 2: q represents the body’s response to this coupling. 
Following this we find a dipolar deformation that results from the coupling between \ a and B a b■ The presence of a 
dipole indicates that the body is accelerated in the tidal field, the acceleration being measured by a a = — MB a bX b > 
this is the Mathisson-Papapetrou spin force [4ll - |43j discussed in detail in Paper 0. And finally, the last set of terms 
describes an octupolar deformation that also arises from the coupling between x a and B a b ; the decaying piece involving 
ft° represents the body’s response. The time-time component of the metric features only the Love numbers ATf 1 , 2i q > 
and .ft 0 ; the remaining Love numbers appear in the remaining components of the metric. 

The light-cone conditions placed on the metric of Sec. IIIII only partially determine the (i/,r, 0, (j>) coordinates, and as 
a result, the metric features a number of arbitrary constants that serve to further specify the choice of gauge. While the 
coordinate system is still geometrically meaningful, the residual gauge freedom implies that the uninteresting gauge 
parameters must be determined alongside the physically interesting Love numbers, making the task more cumbersome 
than it has to be. (An instance of this additional burden can be observed in Ref. Q.) To facilitate the determination 
of the Love numbers in Paper II, in Sec. nv] we present another version of the external metric, this time adopting the 
standard Boyer-Lindquist (i, r, 9 , (f>) coordinates for the background metric, and the familiar Regge-Wheeler gauge for 
the perturbation. 

As emphasized previously, the complete collection of Love numbers {Kf, A"“ as , <£ q , 3°, 25 q , A°) for a selected stellar 
model must be determined by matching an internal metric to the external metric provided in Sec. IIIII or Sec. IIVI For 
a given equation of state for cold matter — the functions p(p) and e(p), in which p is the rest-mass density, p the 
pressure, and e the density of internal (thermodynamic) energy — the Love numbers can be expressed as functions of 
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M/R, the stellar compactness. While it is known that Kf scales as ( R/M ) 5 and A'™ ag scales as ( R/M ) 4 — modulo 
relativistic corrections of order M/R and higher — our considerations in this paper give us no information regarding 
the expected dominant scaling of the remaining Love numbers with R/M. In order to gain some insight into this 
matter, in Sec. |V]we exploit post-Newtonian methods and attempt to calculate the Love numbers of a rigidly rotating 
ball of incompressible fluid. Based on this calculation, we conclude that ~ (R/M) 3 and 5° ~ (R/M) 5 for a generic 
body with an arbitrary equation of state. Our results for the remaining Love numbers are more tentative. Here we 
conclude that <B q ~ (R/M) 3 and £° ~ (R/M) 4 , with an admission that one of these relations might be off by one 
power of R/M; it is thus possible that *B q actually scales as (R/M) 2 , or that A° actually scales as (R/M) 3 . These 
provisional assignments will be confirmed once a proper matching with an internal metric is carried out in Paper II. 

As this work was reaching completion we were made aware of an independent effort [44j to describe the external 
geometry of a slowly rotating material body deformed by tidal forces. This work generalizes ours in the sense that it 
constructs the external metric to second order in the dimensionless spin X a . It is a ls° a restricted version of our own 
efforts, in the sense that it allows only for axisymmetric tidal environments with vanishing B a b- 


II. TIDAL POTENTIALS 

The construction of tidal potentials is presented in great detail in Sec. II of Paper 0 [T}. Here we summarize the 
main results, and introduce the new moments £ a b and B a b that were missed in the earlier work. 

The potentials are obtained by combining Xa, £ ab , B a b , and tt a := x a /r = [sin 6 cos <f>, sin 9 sin (j>, cos 6] in various 

irreducible ways, with each potential carrying a specific multipole order t and a specific parity label (even or odd). 
The coupling of Xa and £ a b produces the odd-parity tensors 

Ra - = £ab\ , Rabc := £(abXc), (2.1) 

in which the angular brackets designate the operation of symmetrization and trace removal, so that R a b c is a symmetric- 
tracefree (STF) tensor. It produces also the even-parity tensor 

Lb := 2 X c e cd{a £ d b) , (2.2) 

in which e a b c is the antisymmetric permutation symbol. The coupling of Xa and B ab produces the even-parity tensors 

Ra ■= B a bX b , Rabc ■= B( abXc )- (2-3) 


and the odd-parity tensor 


Bab := 2 X c e cd{a B d b) . 


(2.4) 


The independent components of the STF tensors £ a b, B a b, Ra, £ab, Rabc, Ra, B a b, and K, a bc can be packaged in 
spherical-harmonic coefficients £ q , £> q , R^,, ££, .F°, £[},, S q , and /C° , respectively. The definitions are summarized 
in Table |U 

The tidal potentials are decomposed in scalar, vector, and tensor spherical harmonics, functions of the angular 
coordinates 9 A = (9, (/). The decomposition involves the scalar harmonics of Table HT1 the even-parity harmonics 
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and the odd-parity harmonics 
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( 2 . 6 ) 


Here £1 A b = diag(l, sin 2 9) is the metric on a unit two-sphere, and D A is the covariant-derivative operator compatible 
with this metric; e A B is the Levi-Civita tensor on the unit two-sphere ( eg$ = sin0), and its index is raised with 
£l AB , the matrix inverse to f 1 A b- It should be noted that the tensorial harmonics are tracefree, in the sense that 

n AB Y%™ = vl ab x% = o. 

The decomposition of the tidal potentials in spherical harmonics is described by 
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TABLE I. Spherical-harmonic coefficients of STF tensors. The relations between IC^, K. a , and are identical to those 
between T a , and x^m- The relations between B^, B a b, and \B ^ are identical to those between £ a b, and x^-m- Finally, 
the relations between KZ,, K a bc, and \Bm are identical to those between J>,, J- a bc, and xTm- 

Xo = X3 = X 
Xic = Xi = 0 
Xls = X2 = 0 

£S = ±(Tn + £22) 

Ti q c = £13 
£ 1 = £23 
£1 = i(fn - £ 22 ) 

£L = £12 

?o = = -2xf 0 q 

Tt=Ti= X £l 
Tl = T2 = x£l 

£q = \ (T11 + T22) = 0 
4 q c = Tl3 = -Xf? s 
T 1% = ^23 = X £ 1c 

4 q c = |(Til - T22) = - 2 X £l 
£1 = £ 12 = 2 X £l 

Jo = |(^ii 3 +^223) = fxT 0 q 

^ c = i(J- 111 +^ 1 22 ) = -^xT 1 q c 
JT s =i(JFn 2 +^ 222 ) = -^xTi q s 

Tie = 1(^113 - J-223) = |xT 2 q c 
n s = >123 = |xT 2 q s 
■Tfc = f (-Tin — 3^122) = 0 
F 2s = i(3Ji ia - F 222 ) = 0 
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It should be noted that tensorial potentials are not defined when £ = 1. The relations between ££ and X £^ displayed 
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TABLE II. Spherical-harmonic functions Y tr ". The functions are real, and they are listed for the relevant modes l = 1 (dipole), 
1 = 2 (quadrupole), and l = 3 (octupole). The abstract index m describes the dependence of these functions on the angle 
for example Y e,2s is proportional to sin 2 <f>. To simplify the expressions we write C := cos 9 and S := sin 9. 

y 1 - 0 = c 

y 1,lc = Scos<t> 

Y Ms = Ssin(f> 

Y 2 ' 0 = 1-3 C 2 


Y 2 ’ic 

= 2SC cosp 
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= 3 S 2 C sin 2(f) 
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= S 3 cos 3(f) 


y 3)3s 

= S 3 sin 3(f) 



in Table HI imply that the tidal potentials associated with £ a b can also be expressed as 

= fr A = -xd*e\, £\ b = ( 2 . 8 ) 

Similarly, we have that 

B\ = - X d^ Al B\ b = - X d*B\ B . (2.9) 

In Paper 0 the metric of a slowly rotating, tidally deformed black hole was presented in terms of the ^-differentiated 
potentials, without recognizing that these are in fact associated with the moments of Eqs. (12.21) and m- 


III. EXTERNAL METRIC OF A TIDALLY DEFORMED BODY — LIGHT-CONE GAUGE 


The external metric of an isolated, slowly rotating body of mass M and dimensionless spin x := |x a | 1 can be 
expressed as 


ds 2 


—f dv 2 + 2 dvdr + r 2 d£l 2 


J\M 2 

r 


sin 2 9 dvd4 >, 


(3.1) 


where / := 1 — 2 M/r and dtt 2 := £lABdO A dd B := dO 2 + sin 2 9 dcf> 2 . The metric is displayed in coordinates (v, r, 9 1 (j>) 
that are tied to the behavior of incoming null geodesics that are tangent to converging light cones. It is easy to show 
that each surface v = constant is a null hypersurface, and that its null generators move with constant values of 9 and 
</>; — r is an affine parameter on each null geodesic. 

The metric of a slowly rotating body immersed in a tidal field produced by remote matter is obtained by perturbing 
the background metric of Eq. Ell). The methods to construct the perturbation are described in detail in Paper 0 (H , 
for the specific case in which the body is a black hole. We adopt these methods here (with very few details provided), 
and make the required changes to account for the presence of matter and the absence of an event horizon. 

In this section we continue to work in light-cone coordinates, so that the coordinates (v, r, 9, cp) keep their geometrical 
meaning in the perturbed spacetime: v continues to be constant on each converging light cone, 9 and <p continue to be 
constant on each generator, and —r continues to be an affine parameter on each generator. These requirements imply 
that g vr = 1, g rr = 0 = g r A , so that g vv , g vr , g v A, and gAB are the only nonvanishing components of the metric. 

The perturbed metric is written as 


gw = — f — r 2 e\ £ q + r 2 e\ £ q + r 2 kf K, d — r 2 k° TC° 


(3.2a) 







6 


TABLE III. Radial functions appearing in the metric of Eq. (13^1) . expressed in terms of x := r/( 2 M), / := 1 — 1/x , and a 
number of integration constants. All functions within square brackets behave as 1 + 0(l/x) when x 1. 

el = f + £[-30x 3 (x- l) 2 ln/+ |x(2x - 1){1 + bx - 6x 2 )\Kf 
e q = / - ^ [2Qx A (x - 1) In / - f x(l + 2x + 6x 2 - 12s 3 )] Kf 
4 = 1 ~ i + Jr [—15 s 3 (2s 2 - 1) In / + 5x 2 (l - 3s - 6 s 2 )] Kf 

b\ = f - Jr [20s 4 (s - 1) In / - §s(l + 2s + 6 s 2 - 12s 3 )] A™ ag 
4 = 1 - i + Jr [—15s 3 (2s 2 - 1) In / + 5s 2 (l - 3s - 6 s 2 )"] A 2 mag 

e ] 1 = ^[60s 5 (s- l) 2 (ln /) 2 + 10s 3 (s - 1)(1 - 9s + 12s 2 ) In/- |s(l - 48s 2 + 108s 3 - 72s 4 )]Kf 
+ & [-30s 3 (s - l ) 2 In / + §s(2s - 1)(1 + 6 s - 6 s 2 )] & + f 2 & - ^ 7 q + ^ ~ ^ + si* 

[-^s 6 (s - l)(ln /) 2 + f s 4 (l + 10s - 12s 2 ) In/ + f s 2 (2 + 9s + 24s 2 - 36s 3 )]Kf 

- £; [ 20 s 4 (s - 1 ) In / - §s(l + 2 s + 6 s 2 - 12 s 3 )] & + §/£ q + ^( 2 s + l) 7 q + ^ 

e q 7 = [20s 5 (2s 2 - l)(ln /) 2 - f s 4 (l + 4s)(1 - 3s) In/ + ^s 2 (l - 6 s 2 + 36s 3 )]Kf 

+ sfs [~15x 3 (2x 2 ~ 1) In / + 5s 2 (1 - 3s - 6 s 2 )] £ q + ±(1 - ^)e q + ^ 7 q + ^ 

b\ = [-M2s 6 (s - l)(ln /) 2 + f s 4 (l + 10s - 12s 2 ) In / + fx 2 (2 + 9s + 24s 2 - 36s 3 )] A'“ ag 

- i [ 20 * 4 (* - 1 ) In / - §s(l + 2 s + 6 s 2 - 12 s 3 )] © q + f /< 8 q + ^- ^- 3 ^-^ 

6 ? = - J i 7 [20s 5 (2s 2 - l)(ln /) 2 + f s 4 (4s + l)(3s - 1) In / + f s 2 (l - 6 s 2 + 36s 3 )] K™ e 

+ 35* [~15x 3 (2x 2 - 1) In/ + 5 s 2 (1 - 3s - 6 * 2 )]< 8 q + ±(1 - ^)® q + 4 ^ - 

ki = -g|r [15 s 3 (s - l)(2s - l)(5s - 1) In/ - fs(l + 2s - 34s 2 + 144s 3 - 120s 4 )] A“ ag 
+ Xc d + L- 4 —4_ 

16a? 4 x 10a? 2 5o? a 

fc 4 = sYT [-30S 5 (5s - 4) In / - |s 2 (2s - 1)(1 + 6 s + 30s 2 )] A'” ag 

1 / I J_ _ 2 _ c 

16a: 4 4” 2a; 5a? 20a: 4 

k\ = 4j. [—10s 3 (s - 1)(1 + 3s + 140s 2 - 420s 3 + 280s 4 ) In / + fs(l + 12s + 34s 2 + 244s 3 - 3640s 4 + 6720s 5 - 3360s 6 )] A'” ag 

+ ^ [—420s 4 (2s - 1)(s - l ) 2 In/ + 7s 2 (l + 10s - 130s 2 + 240s 3 - 120s 4 )]H° + 2/ 2 (l - ±)x$° + ^c 0 + ^ ^ 

k% = iffr [|§s 4 (5 — 4s + 280s 2 - 700s 3 + 420s 4 ) In/+ §s 2 (7 + 32s + 46s 2 + 420s 3 - 2940s 4 + 2520s 5 )] A“ ag 

+ Jr [210 s 5 (3s - 2)(s - 1) In/ + |s 2 (l + 5s + 30s 2 - 210s 3 + 180s 4 )]- f/(I - £)x&° + §§^c° + ^ 

fe° = - [-f s 4 (l + 12s - 140s 3 + 140s 4 ) In/ + f s 2 (l - 2s 2 + 140s 3 + 420s 4 - 840s 5 )] A'” ag 

- £ [—84s 4 (1 - 10s 2 + 10s 3 ) In/ - 14s 3 (l - 10s - 30s 2 + 60s 3 )]- (1 - ± + t ^ T )s^° + j^c° + 

ft = - s yr [—30 s 5 (5s - 4) In/ - §s 2 (2s - 1)(1 + 6 s + 30s 2 )] A 2 el + + ^ 7 d + ^ 

ft = -^[|s 3 (4s - 5) In / - |s(7+ 18s - 24x 2 )]Kf 

+ [210s 5 (3s - 2)(s - 1) In/ + |s 2 (l + 5s + 30s 2 - 210s 3 + 180s 4 )]5° + |/(1 - £)xd° + 5 ^ - 

f? = - [- 6 s 3 (2s - 1 ) In / - s(l + 12s 2 )] Kf 

+ £ [—84s 4 (1 - 10s 2 + 10s 3 ) In/ - 14s 3 (l - 10s - 30s 2 + 60s 3 )]S° - (1 - | ^)*5° + ^ 7 ° + ^ 


g vr — 1, 


9 M 2 9 A , 

9vA = —^tXa - 3 ( e 4 £a - b 4 B a) + 4 {$Pa - K B a) - 4 

g AB = r 2 n AB - ir 4 (e? - 6? b< AB ') + 4 (e? - 6? 


(/ 4 d - fc d /Cl) + r 3 (/| 71 + fc 4 ° , 

-r\f° 7 J° AB -k° 7 1C AB ), 


(3.2b) 
(3.2c) 
(3.2d) 


in which e q , b ^ e^, 6^, fc d , fc°, / d , and f° are functions of r that are determined by solving the vacuum Einstein field 
equations. The radial functions are listed in Table uni 

The strategy to solve the held equations is as follows. First, we switch perspectives and consider the metric of 
Eq. (13.21) to be a perturbation of the Schwarzschild solution expressed in (v, r, 9, <f>) coordinates — this is Eq. (13.11) 
with x se t equal to zero. 

Second, we manufacture a first-order perturbation to this new background metric by introducing the rotational 
potential x A an( t the purely tidal potentials constructed from E a b and B a b- With the spherical harmonic decompositions 
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of Eqs. <H3J|, the linearized field equations reduce to a system of coupled differential equations for the radial functions 
inserted in front of these potentials. Schematically, these take the form of 

^>1 = 0, (3.3) 

where Jzf'j. is a second-order differential operator, and wf(r) is the collection of radial functions that appear in the 
first-order metric perturbation. (In practice, the equations for the rotational, £ = 1 perturbation decouple from 
the equations for the tidal, l = 2 perturbations, and the equations for {e q ,e q ,e q } decouple from the equations 
for { 64 , fry}.) Equation (13.31) is then integrated, and the general solution is formed from two linearly independent 
modes, one decaying as r increases, the other growing. In the case of the rotational perturbation, the decaying mode 
produces the angular-momentum term in Eq. (EH), and the growing mode is set to zero, because it corresponds to 
an uninteresting coordinate transformation to a rotating frame. In the case of the tidal perturbation created by £ a b 
(respectively B a b ), the growing mode represents the external tidal field, and the decaying mode represents the body’s 
response to this tidal field, measured by the gravitational Love number Kf (respectively AT™ as ), which appears in the 
solution as an integration constant. The solution for b q is also observed to involve c, an additional integration constant 
that represents a residual gauge freedom that preserves the geometrical meaning of the (v, r, 0 , <f>) coordinates. 

Third, the first-order perturbation is used as a seed to construct a second-order perturbation that accounts for the 
coupling between x°' and Sab , and the coupling between \ a and B a b . This perturbation ignores terms quadratic in x“, 
and terms quadratic in £ a b and B a b . The composition of the l = 1 and 1 = 2 spherical harmonics produces the dipole, 
quadrupole, and octupole potentials listed in Sec. El and insertion of the radial functions gives rise to the metric of 
Eq. (13.21) . Substitution of this metric in the vacuum field equations yields differential equations of the schematic form 

^> 2 fe = (3.4) 

where w 2 is the collection of radial functions that appear in the second-order perturbation, S J is a set of source terms 
constructed from the first-order radial functions, and J£ >j k is the same differential operator as in Eq. m ■ The general 
solution to Eq. EH, 


w 2 = ^(particular) + ^(decaying) + ^(growing), (3-5) 

is a linear superposition of a particular solution to the system of differential equations, a decaying solution to the 
homogeneous system C\w\ = 0, and a growing solution to the same homogeneous systenfl 

In Table Ell the decaying solutions to the homogeneous system are identified with the integration constants <£ q , 
3°, $ 8 q , and M°, and the growing solutions are identified with £ q , 3 d ’°, ® q , and Jt°. In addition to these, the solutions 
depend on gauge constants 7 d,q, ° and c d,q, ° that (like c) specify the residual freedom of the light-cone gauge. It is 
useful to note that with Kf = A“ ag = (£ q = 3 ° = <B q = 8° = & = 3 d ’° = % q = £° = 0, the solutions of Table Hill 
reduce to those obtained in Paper 0 Jl|. 

Many of the integration constants introduced in Table Mil do not have a physical meaning, and can be set equal to 
zero without loss of generality. We have already put in this category the gauge constants c, 7 d,q ’°, and c d,q ’°, which 
merely serve to further specify the coordinate system. We can also single out as unphysical the constants £ q , 3 d ’°, 
S' 1 , and j| 0 , which come with the growing modes of the radial functions. 

Let us first examine the constant 3 d , which appears in / d . It is easy to see that this term can be removed from g v A 
by making the replacement 

X a - 8M 2 5 d T“ -> X a (3-6) 

in the metric. This shift of the original spin vector is unobservable, and its effect on other terms in the metric scale 
as x 2 and can be neglected. Thus, we see that 3 d can be set equal to zero without producing a different physical 
situation. 

Turning next to the terms involving the constants C: q and *B q , we see they also can be removed by performing the 
shifts 


£ab - £ q £ab £ab , & - & Kf -> C: q 


(3.7) 


3 The growing and decaying solutions can be unambiguously identified. The decaying solution decreases with increasing r, faster than any 
term in the growing solution, and it is singular in the limit r —> 2 M. The growing solution increases with r, and is regular in this limit. 
It should be pointed out that the limit r —>■ 2 M is entirely mathematical, since the presence of a meterial body inside r = R > 2 M 
implies that r = 2 M lies beyond the metric’s range of validity. The mathematical limit is nevertheless available to identify the growing 
and decaying solutions unambiguously, and the fact that it cannot be realized physically is of no concern to this identification. 



and 


B ab -WB ab ^B ab , 25 q — ® q /\™ as —t ® q - (3.8) 

These shifts also are unobservable, and their effects on other terms in the metric can also be neglected. We therefore 
conclude that £ q and <B q can be set equal to zero without loss of generality. 

The case of the constants ft° and 3° is more subtle, because the metric of Eq. (EH does not include a tidal 
perturbation generated by octupole moments £ abc and B ab c that could be shifted by rotational corrections. But this 
situation merely reflects an incompleteness of our description, which can easily be supplemented with the missing 
octupole terms (see Refs. (HI, Hi}). With the octupole contribution included, we find that the terms involving j|° and 
3 ° can be eliminated by performing the shifts 

M£ abc + 3 8°lC abc -> M£ abc , (3.9) 

and 

\MB abc + 3Z°T abc -> ^MB abc , 3° + 23°AT ag -a 3°, (3.10) 

where I\f and A'“ ag are the £ = 3 gravitational Love numbers. The factor of | in front of B abc reflects the definition 
of the £ = 3 tidal potentials adopted by Poisson and Vlasov [if| ; its original source is the choice of normalization for 
the tidal moments made in Ref. [47|. Once again the conclusion is that R° and 3° can be set equal to zero without 
altering the physical situation. 

With the constants c, 7 d,q '°, c d,q, °, £ q , 3 d ’°, ® q , and A° dismissed as physically uninteresting, we are left with 
a metric that depends on the two original Love numbers ATf 1 and AT™ ag , as well as the remaining four integration 
constants <£ q , 3°, ® q , and A°. These are attached to the decaying solutions of the homogeneous perturbation equations, 
and we interpret them as a new class of Love numbers associated with the coupling between the body’s rotation and 
the external tidal field. Fittingly, we may refer to l£ q , 3°, S q , and J?° as rotational-tidal Love numbers. 

The new Love numbers are thus identified as integration constants in front of the decaying solutions to the homo¬ 
geneous perturbation equations. They are gauge-invariant: While the precise expression of the decaying solution will 
be altered by a change of gauge, its identity as a decaying solution will be preserved, and the numerical value of € q , 
3°, ® q , and A° will also be preserved by the gauge transformation. The gauge-invariant nature of the Love numbers 
can also be recognized from a computation of gauge-invariant quantities. For example, the metric of Eq. EH) can be 
used to calculate the Newman-Penrose quantity ip o in a null tetrad adapted to the principal null congruence of the 
slowly rotating background spacetime; as is well known, ipo vanishes in the background spacetime and is therefore 
gauge-invariant in the perturbed spacetime. We carried out this computation but choose not to disclose the results 
here to avoid displaying the long expressions. We state nevertheless that ipo is observed to depend on £ q , 3°, 23 q , 
and A°, which each constant multiplying a linearly independent solution of the Teukolsky equation. Because each 
term is linearly independent and therefore individually gauge-invariant, we can conclude that the rotational-tidal Love 
numbers are indeed gauge-invariant. 

Why four new Love numbers? As we have seen, the counting of Love numbers is governed by the coupling between 
spherical harmonics when the £ = 1 rotational perturbation is combined with the £ = 2 tidal perturbation in a 
second-order computation of the complete perturbation. The coupling produces terms of multipole orders £ = 1, 
£ = 2, and t = 3, and each term comes with a number of essential integration constants, which are distinguished from 
those associated with the residual gauge freedom. In the case of £ = 2, the number of essential constants is two for 
the even-parity sector, plus two for the odd-parity sector, and we have seen that two of these can be absorbed into a 
redefinition of the tidal quadrupole moments; the two remaining constants are the £ = 2 rotational-tidal Love numbers. 
The conclusion is the same for £ = 3, and we have our total of four new Love numbers. The dipole case is exceptional, 
and we have seen that here, the single essential constant can be absorbed into a shift of the angular-momentum vector. 

We may generalize the argument and conclude that the coupling between the £ = 1 rotational perturbation and a 
tidal perturbation of multipole order £ > 3 will produce six new Love numbers. In this case the coupling produces 
terms of multipole order £ — 1 , £, and £ + 1 in both the even-parity and odd-parity sectors, with each multipole of 
each sector coming with two essential integration constants. One of these can be absorbed into a redifinition of the 
associated tidal multipole moment, and we are left with 3 multipoles x 2 sectors = 6 new Love numbers. 


IV. EXTERNAL METRIC OF A TIDALLY DEFORMED BODY — REGGE-WHEELER GAUGE 

In this section we repeat the calculations carried out in the preceding section, but express the metric perturbation 
in the Regge-Wheeler gauge instead of the liglit-cone gauge. And instead of the light-cone coordinates (u,r, 9, 0lc)> 
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TABLE IV. Radial functions appearing in the metric of Eq. 114.411 . expressed in terms of x := r/(2M), f := 1 — 1/x, and a 
number of integration constants. All functions within square brackets behave as 1 + 0(l/x) when 1. 

4 = f + Jr [—30a: 3 (a; - l) 2 In / - §*(2® - l)(6x 2 - 6a: - l)]AT 2 el 

p q _ f- 2„q 

e q = 1 - 2 ^ + pr [-15a: 3 (2a: 2 - 1) In / - 5a: 2 (6a: 2 + 3a; - 1)] AT 2 el 
= / - J. [20x 4 (x - 1) In / + §x(12x 3 - 6x 2 - 2a; - 1)] AT 2 mag 


4 = Jr [—30a: 3 (a; - l ) 2 In / - §x(2x - l)( 6 x 2 - 6x - l)](£ q - j^g) + f 2 £ q 
4 = - 3 ^ [15a: 4 (3a: - 1) In / + f^bx 3 - 30x 2 - 1)] K? + & - ^ 

sq _ 2 -q 

'-'7T J '-'tt 

e q = Jr [—15x 3 (2x 2 - 1 ) In / - 5x 2 (6x 2 + 3x - 1 )] (<f q - + 5 ^ ( 2 a : 2 - l)£ q 


6 q = - Jr [20x 4 (x - 1) In / + §x(12x 3 - 6x 2 - 2x - 1)] (23 q - ^ + §/© q 

Lq _ 3 r_ 20 X 5 (2x — 1) 1 r _ 10 x 2 (12x 3 +x+l) 1 rv-mag , 1 

r Ax 7 L 3 x —1 J 9 x —1 J XV 2 ' IS 


12 ;r 2 (x —1) 


1 20x 2 -49x+38 

20 x 2 (x —1) 


4 = -s|t [f (20x - 9)(x - l)x 4 In/ + |-^(240x 5 - 468x 4 + 242x 3 - 16x 2 - x + 2)]V 2 mag + T^bryC d + 

4 = ■£? [f 1 n/ + |(^ F (12 x 4 - 18x 3 + 4x 2 + x + 2)]^ + 
fc? t = Jr [—10x 4 (x - l)(280x 3 - 420x 2 + 140x + 3) In / - 2800x 7 + 5600x 6 - + 2I°x 4 + ±fx 3 + 5x 2 - |x - f] AT? ag 

+ Jr [— 420x 4 (2x - l)(x - l) 2 In/ - 7x 2 (120x 4 - 240x 3 + 130x 2 - lOx - 1)]£° + / 2 (2x - 1)1° + ^ ^ 

k° rr = lr [-12 (280x 3 - 420x 2 + 140x + 1) In / - J ^4^ (3360x 7 - 6720x® + 3640x 5 - 268x 4 - 34x 3 - 2x 2 + 3x - 5)] A'” ag 

+ Jr [— 420x 6 (2x - 1) In / - 7 I ^ [F (120x 4 - 240x 3 + 130x 2 - lOx - 1)] £° + (2x - 1)6° + 
k° = [—10x 5 (140x 3 - 140x 2 + 13) In / - |x(1680x 6 - 840x 5 - 280x 4 + 16x 3 - 6x 2 - 4x - 1)] A'” ag 

- Jr [—84x 4 (10x 3 - 10x 2 + 1) In/ - 14x 3 (60x 3 - 30x 2 - lOx + 1)]£° - gjr(10x 3 - 10x 2 + 1)R° + gj- 


ft = -28^ [—84 x 4 (10 x 3 - 10x 2 + 1) In / - 14x 3 (60x 3 - 30x 2 - lOx + l)]K? + £$ d + ^+2^ 
f? = -^[-|a: 2 (5x - 2) In / - |x(10x + 1)] K? 

+ Jr [210x 5 (3x - 2)(x - 1) In/ + |x 2 (180x 4 - 210x 3 + 30x 2 + 5x + 1)]3° + ±/(3x - 2)5° - ^ + gj- 


we cast the background metric in the standard Boyer-Lindquist coordinates (t, r, 9 , </>), so that the line element is now 
expressed as 

ds 2 = —/ dt 2 + f^ 1 dr 2 + r 2 d£i 2 — 2— -sin 2 9 dtd(/>, (4-1) 

r 

where / := 1 — 2 M/r and d£l 2 := flABd9 A d9 B := d9 2 + sin 2 9 dip 2 . It is important to note that the coordinate 4> 
adopted here is distinct from the </lc featured in Sec. IIIII (and denoted cj) there); we have the relations dv = dt + f~ 1 dr 
and define = d<j>+ (2x'M 2 /r 3 /) dr between the coordinate systems. 

The Regge-Wheeler gauge conditions are formulated at the level of each £m mode of the metric perturbation. For 
the even-parity sector generated by £ a b , £ab, /C a , and K, a bc , we decompose the metric perturbation p a p as 

Pab = J2 h *™ Y " m ’ P«B = Y J i t ™ Y A, PAB=r 2 J2{K em ^ABY £m + G em Y^), (4.2) 

£m £rr\ 

where x a = ( t,r ) and 9 A = (9,<j>), and we impose 4 m = 0 = G im for i = 2 and £ = 3. For £ = 1. the tensorial 
harmonics Yfy and the associated perturbation variables G lm are not defined, and we adopt instead A' 1,m = 0 as a 
gauge condition, along with j 4 ’ m = 0. For the odd-parity sector generated by B a b, Bab , -4, and £F a bc, we decompose 
the metric perturbation as 

Pab = 0, PaB = J2 h * mX A, PAB = Y J ht 2 mX AB (4-3) 

irr\ Im 

and impose = 0 for £ = 2 and £ = 3. For £ = 1, and hf™ are not defined, and we adopt instead h 4 ’ m = 0 as 
a gauge condition; this choice is motivated by the fact that h 3,rn is found to vanish by virtue of the field equations. 
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The strategy to integrate the vacuum Einstein field equations is readily adapted from Sec. IIII1 and we obtain a 


perturbed metric of the form 

9tt = -/ - r 2 4 £' 1 + r 2 e q tt £ q + r 2 k d tt lC d - r 2 k° tt 1C°, (4.4a) 

g tr = r 2 e q tr £ q 1 (4.4b) 

9rr = /- 1 - r 2 e q r £ q + r 2 e q r £ q + r 2 k d r /C d - r 2 k° r /C°, (4.4c) 

gtA = -^-X d A + f r 3 6j 6^ - r 3 b q - r 3 f d T d A + r 3 f? T° A , (4.4d) 

grA = —r 3 b q B\, (4.4e) 

g AB = r 2 n AB - r 4 e q £% + r 4 e q £l AB £ q + r 4 k° V AB 1C°, (4.4f) 

with the radial functions listed in Table hyi 


The radial functions depend on a number of integration constants, but many do not have a physical meaning. As 
was discussed in Sec. m the terms involving £ q , 5 d '°, * 8 q , and A° can be eliminated by redefining x“> £ a b, B a b, 
£ a bc, and Babe] these constants can therefore be set equal to zero without loss of generality. The constants c d and y d 
specify the residual gauge freedom associated with the 1=1 even-parity and odd-parity modes, respectively. (There 
is no residual gauge freedom in the t > 2 inodes; in this case the Regge-Wheeler conditions completely specify the 
coordinate system.) 

The remaining integration constants,£ q , 5°, ® q , and A° can again be interpreted as a new class of rotational-tidal 
Love numbers. Computation of ip o from the metric of Eq. (14.41) and comparison with the results obtained (but not 
displayed) in Sec. IIIII confirms that the Love numbers introduced here are the same as those introduced in Sec. m 
This consistency check explains the curious shifts £ q —> (£ q — ^, QS q —!► *B q — o bserved in Table [IVJ without these 
the metric of Eq. (14.41) would not be directly obtained from the metric of Eq. (13.21) after a coordinate transformation 
followed by a gauge transformation. 


V. SLOWLY ROTATING, TIDALLY DEFORMED BODY IN POST-NEWTONIAN THEORY 

In this section we develop a post-Newtonian description of the geometry of a slowly rotating, tidally deformed 
material body. We perform a leading-order calculation that accounts for the coupling between \ a an d £ab ; the 
coupling between x“ and B a b occurs at a higher post-Newtonian order and is neglected here. We adopt the simplest 
model for the slowly rotating body: an incompressible fluid with uniform density po , rotating rigidly with a constant 
angular velocity oj a . Our main objective in this section is to obtain useful information regarding the scaling of (£ q , 
5°, 23 q , and A° with the body’s radius R. 


A. Newtonian structure 

The unperturbed body has a density function p = poQ(R — r), with R denoting the body’s radius, an internal 
potential U- m = (GM/2R)(3 — r 2 /R 2 ), an external potential U ou t = GM/r , and its pressure profile is given by 
p = (3GM 2 /8 ttR 4 )(1 — r 2 /R 2 ). The body’s mass is related to its density by M = |7rp 0 R 3 - 

The tidal deformation produces perturbations Sp , SUi n , 5U ou t, and Sp to these quantities, as well as a perturbation 
SR to the position of the boundary. With a tidal environment characterized by a quadrupole moment £ a b , we have 
that (5f7 ou t can be expressed quite generally as 

SU out = ~[l + 2k 2 (R/r) 5 ]£ ab x a x b , (5.1) 

in terms of a gravitational Love number k 2 ] the first term (which grows as r 2 ) represents the external tidal field, while 
the second term (decaying as r -3 ) represents the body’s response to the tidal forces. With this information, SR is 
deduced from the requirement that the boundary must be an equipotential surface; we obtain 

S R = ~^^(l + 2k2)£ ab n a n b . (5.2) 

The perturbation of the internal potential must have a pure quadrupolar form proportional to £ a bX a x b , and the 
constant of proportionality is obtained by continuity of SU- m and 5U ou t at r = R; we find 

SU- m = -^(1 + 2k 2 )£ a bX a x b . 


(5.3) 
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The density perturbation is obtained from the new assignment p + Sp = po@(R + SR — r). Expansion of the step 
function to first order in SR and use of Eq. (E1D produces 


Sp=~ — (l+2k 2 )£ ab n a n b S(r-R). (5.4) 

07rG 

To obtain the pressure perturbation we appeal to the statement of hydrostatic equilibrium, d a Sp = pd a SU ln + 8pd a U m . 
The angular components of this equation integrate to Sp = poSU- ln , and we get 

Sp = - (1 + 2k 2 )£ ab x a x b . (5.5) 

It can be verified that p + Sp properly vanishes when r = R + SR. The results for SR, Sp, and Sp can alternatively be 
derived on the basis of the Lagrangian displacement vector 

R 3 

^ a = ~2 Gl (1 + 2fc2) ^' (5-6) 

Finally, k 2 can be determined from the discontinuity in d r 5U at r = R implied by Poisson’s equation V 2 5?7 = — AnGSp 
and the 5-function of Eq. (15.411 . A short calculation reveals that the discontinuity must be given by [d r 5U] = 
|(1 + 2k 2 )R£ ab £l a £l b , where [d r SU] := d r U ou t — d r U m evaluated at r = R. Substituting our previous results for 5U m 
and SU ou t, we arrive at 

k2 = \, 1 + 2*2 = ! (5J) 

This is the well-known value of the gravitational Love number for an incompressible fluid. 

The body is taken to be rotating rigidly with a uniform angular velocity ui a . The rotation is assumed to be slow, 
and we work consistently to first order in uA The unperturbed velocity field inside the body is v a = e a bcOJ b x c , 
the spin angular momentum of the unperturbed body is given by S a = jrMR 2 oj a , and its dimensionless version is 
X a = ( c/GM 2 )S a , or 


In terms of this the velocity field becomes 


X 


2 cR 2 
5 GM 


to . 


V a 


5 GM b c 

2^ eabcX X 


(5.8) 


(5.9) 


The tidal deformation creates a perturbation Sv a = v b d b £, a — £, b d b v a of the velocity field, which can be calculated from 
the Lagrangian displacement vector of Eq. m - We obtain 

pc. 73 „ 

Sv a = - £ ab X b . (5.10) 

8 c 

The tidal perturbation changes the body’s angular momentum by 5S a = e ab c f x b 5j c d 3 x, where Sj a = Spv a + p5v a - 
A short calculation reveals that the shift in the dimensionless spin is given by 


fix' 


5 R 3 


AGM 


£ a b x b 


5 R 3 


AGM 


T a 


(5.11) 


B. Vector potential 


Our post-Newtonian calculation is based on the metric 


goo — ~~ 1 H — 2 U + 0(c ), gQ a — — -U a + 0(c ), g ab — 


1 + ~^ u 

& 


+ 0(c 4 ), 


(5.12) 


which involves a vector potential U a in addition to the Newtonian potential U encountered previously. The coupling 
between x a smd £ ab is captured by the vector potential, which satisfies the field equation 


V 2 1/ q = -AnGpv a . 


(5.13) 
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An examination of the post-Newtonian equations (see, for example, Sec. 8.1 of Ref. [48|) reveals that the coupling does 
not appear in g 00 at first post-Newtonian (Ipn) order, and we therefore omit the Ipn terms at order c~ 4 . Similarly, 
consideration of g a b at 2pn order reveals that the coupling does not appear within the omitted c~ 4 terms. The 
unperturbed vector potential is 


outside the body, and 


U° 


{GM) 2 
2 cr 3 


C abcX X 


U" 


{GM) 2 
4 cR 3 


(5 - 3 r 2 /R 2 )e abc x b x c 


(5.14) 


(5.15) 


inside the body. The additional factor of c -1 results from expressing the angular-momentum vector in terms of its 
dimensionless version of Eq. (15.81) : this promotes go a to a quantity of 1.5pn order. 

The perturbation of the vector potential is sourced by the sum of 6\j a := Spv a and S 2 j a = pSv a , and correspondingly 
it is expressed as the sum of 


SiU a (x) = gJ 


S da(x') 3 r tt - n [ 5 ^a{x') ,3 , 


\X ~ X' 

To evaluate this we rely on the addition theorem 


5 2 U a {x) = gJ- 


x 


^ = E 2^1 

tm > 


(5.16) 


(5.17) 


the identity [Eq. (1.171) of Ref. [48j; L := aia 2 ■ ■ ■ ai is a multi-index that includes a number £ of individual indices, 
Q l := fl ai Q a2 ■ ■ ■ fl ae , and the angular brackets indicate that all traces are to be removed from fl L ] 


f Yl m {0',<j>') QfW dfl' = (5.18) 


and a decomposition of 

:= e abc x b £d e n c n d n e (5.19) 

into irreducible components. This is accomplished with the identity l e = + ^{8 cd VL e + 5 ce Q d + S de fl c ), 

which gives rise to Z a = |Z 4 + Zf, where 

:= e abcX b £ c d n d , Zl := e abcX b £ de & cde l (5.20) 

These quantities can be expressed in terms of the tidal potentials [l| 

Ti := \e abc £l b T c , T° a := \e abc £l b T c de £l d £l e , (5.21a) 

£ q := ^£ ab n a n b , £?, := 1 ( S b - fl a n b ) £ bc £l c , (5.21b) 

which are Cartesian versions of the angular potentials introduced in Sec. El the relation is given, for example, by 
T\ = where fl a A := dQ a /dd A . We have 

c- 2 z a 4 = \{K + £2 + n a £') , c~ 2 zl = T (2£2 - 3f2 a f q ) - r a . ( 5 . 22 ) 

The computation of <5i U a and S 2 U a proceeds by inserting 

<5i ja = ~-^- — Z a S{r - R), S 2 j a = ^-^£ ab x b Q{R- r), (5.23) 

32-7T C 327T cR 

as well as the addition theorem of Eq. (15.171) , within Eq. (15.161) . The angular integrals are evaluated with the help of 
Eq. (15.181) . after involving the decomposition of Z a into Z 4 and Z^. The final expressions, 


5iU° a ut 


-5 GMc 


' 1 R 3 

15 'r 2 



1 R 5 

= ^£ 3 -^l n £ 3 -T° 


7 r 4 V 15 


(5.24) 
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and 


71 

S 2 u° ut = —GMc 




are obtained after making use of Eq. (15.2211 . The complete perturbation is 


75 

6U° ut = —GMc 


1 R 3 


1 R 5 


1 R 5 


15 ^ + 77^( 2 4 q - 30^) ^ 


105 r 4 


7 r 4 


(5.25) 


(5.26) 


and we note that the terms proportional to ( R 3 /r 2 )£ a b that appear in both 6iU° ut and S 2 U° ut cancel out after taking 
the sum. We note also that the term proportional to T 3 can be fully accounted for by the shift in spin vector 
described by Eq. (15.111) . With the factors of c -2 contained in the tidal potentials, we find that SU° ut scales as GM/c 
and contributes to g Qa at 1.5 pn order. 


C. Light-cone coordinates 

The post-Newtonian metric of Eq. (15.121) cannot be compared directly with the metrics obtained in Secs. m and 
hyi because they are expressed in different coordinate systems. (The potentials U and U a are now the full, perturbed 
potentials; they were previously denoted U + SU and U a + 5U a , respectively.) Here we transform the post-Newtonian 
metric to the light-cone coordinates of Sec. 1III1 and compare the results with the metric of Eq. (13.21) . For this purpose 
it is useful to look at the post-Newtonian metric as a perturbation of the Minkowski metric, write g a p = rj a p + h a p, 
and work consistently to first order in the perturbation h a p. 

The first step is to perform the coordinate transformation 

ct = cv — r, x a =rn a (d A ), (5.27) 

which implies that v is a null coordinate in flat spacetime. The transformation brings the Minkowski metric to the 
form 


Voo = -1, Vor = 1, Vab = r 2 n A B , 

(5.28) 

and the perturbation to the form 


2 

hoo = — U, 
c- 

(5.29a) 

2 4 

h 0 r = u a n a , 

c c 3 

(5.29b) 

4 

hoA = mG a £l A , 

e 3 

(5.29c) 

4 8 

h rr = -u + -u a n a , 

C" c 3 

(5.29d) 

h r A = a £l Al 

C 3 

(5.29e) 

2 , 

h-AB = ~n r UU, AB . 

<r 

(5.29f) 

The second step is to perform a gauge transformation, 


^ck/ 3 ^ ^a/3 — ^a/3 ^a£/3 

(5.30) 

to ensure that v remains null in the perturbed spacetime. This requires h' 0r = h' rr 
computation reveals that the gauge vector £ a is determined by 

= h' rA = 0, and a simple 

drCo = ~u-^u a n a , 

c z c 3 

(5.31a) 

drtr = ^U+^U a n a , 

(5.31b) 

r 2 d r (r~ 2 £ A ) = \rU a £l A - d A £ r . 

(5.31c) 
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With the null conditions satisfied, the remaining components of the metric perturbation become 

*■00 = ^2 (5.32a) 

h' QA = -^rU a fl A - d A (o, (5.32b) 

h' AB = ^ 1 " 2 UCI A b — D a £ b - D b ^ A - 2r(£ 0 + ^ r )^AB- (5.32c) 

We focus our attention on the piece of the metric associated with SU a , which captures the coupling between \ a and 
Sab- We insert Eq. (15.261) within Eq. (15.311) and obtain 


SS,o 

SSr 

*U 


75 GM ( 

YYY 

75 GM 

YYY 

75 GM 


1 R 5 


V105 r 3 
1 R 5 
105 T 3 " 


£ q +Po 


£ q +P , 


—i? 3 J^ - — F a + rd A p r + r 2 /3 A 


(5.33a) 

(5.33b) 

(5.33c) 


where /3o, f3 r , and f3 A are arbitrary functions of d A . For our purposes here it is sufficient to restrict the residual gauge 
freedom to functions of the form 


Po =P q 0 R 2 S\ p r =p q R 2 £ q , p A =p d RT d A + p°RI 1, (5.34) 

in which p q , p q , p d , and p° are arbitrary dimensionless coefficients. Inserting the gauge vector within Eq. (15.321) 
produces 


bh' A B = 


75 GM 

YYY 

75 GM 


1 R 3 
15 r 


i-—T A A -2piRn\-Y-YT° A 


1 R 5 
7 r 3 


2(2 p q - p q )R 2 rQ AB S q - 2 p q R 2 rS q AB + 


(J_R*_ 

\14 r 2 


2 p°Rr 2 



(5.35a) 

(5.35b) 


We notice that terms proportional to (R 5 /r 3 )E\, which used to appear in Sh 0A , no longer appear in Sh' 0A . We can 
further eliminate the term proportional to f l AB S q in Sh' AB by setting 2 p q = p q ; this constitutes a refinement of the 
light-cone gauge. 

By comparing Sh' 0A and Sh' AB with the metric of Eq. (13.21) . we can read off the post-Newtonian expressions for the 
relevant radial functions. We have 
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Comparing these with the expressions listed in Table cm we find a precise match at this leading, 1.5PN order. Our 
expression for e q matches the j^x _3 7 q term in Table cm and we see that the constant p q can be related to y q . 
Similarly, our e q corresponds to the -^x~ 3r ) q term in Table [Till Our expression for / d , which reflects the shift of the 
body’s angular-momentum vector created by the tidal perturbation, matches the expected £~ 4 5 d term from Table ITTT1 
and we recall that the parameter 5 d was indeed related to a shift of \ a in Sec. Mil Our expression for matches the 
expected 2cc -6 3° term from Table cm and in this case we can assign the precise value 
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to the associated rotational-tidal Love number. This assignment is confirmed by comparing f° with the expression of 
Table Hill and in addition, we see that the gauge constant p° can be related to 7 0 . 

Our calculation cannot produce an expression for e 9 , which occurs at a higher post-Newtonian order, and it also 
cannot produce the terms that couple \ a to B a b in the metric of Eq. (13.21) . In addition, our calculation did not 
produce the expected — 2x~ 5 S q term in e q , nor the expected |a: _5 2 q term in e 9 . The reason for this must be that 
these terms occur at a higher post-Newtonian order. A match at 1.5PN order would have implied a scaling of the form 
(GM/c 1 2 ) 5 £ q oc ( GM/c 2 )R 4 for 0: q . In the absence of such a match, we may expect that the term appears instead at 
2.5PN order, with an expected scaling of (GM/c 2 ) 3 l£ q oc ( GM/c 2 ) 2 R 3 for £ q . 

As a final note, we mention that our post-Newtonian calculation did not produce the expected i;r _1 7 d term in / d . 
The fault here does not lie with a failure to go to a sufficiently high post-Newtonian order. It lies instead with the 
omission of a dipolar solution to V 2 [/ a = 0 as an additional piece of the vector potential of Eq. (15.261) . It is easy to 
show that the homogeneous term 5\JJ a = \t d GMcrR 3 is indeed the required missing piece. 


D. Conclusion: Scaling of rotational-tidal Love numbers 


The consistency between the post-Newtonian metric obtained in this section and the exact metric obtained in 
Sec. m is pleasing and reassuring, but our main purpose was not to demonstrate this consistency. Our goal was 
instead to seek guidance in the scaling of (£ q , 55 q , 5°, R°, the new class of rotational-tidal Love numbers, with the 
body’s radius R. The calculations presented here allow us to anticipate that for strongly self-gravitating bodies with 
arbitrary internal structure, 5° will definitely scale as R 5 , and £ q will likely scale as R 3 . 

Our calculations give us no direct guidance regarding A° and Q3 q , but inspection of the equations of post-Newtonian 
theory (including the statement of hydrostatic equilibrium) indicates that at leading order, the coupling between \ a 
and B a b produces terms of order c~‘ in both goo and g a b ■ Matching this observation with the x~ 6 A° terms in Table Hill 
produces the expectation that M° should scale as (GM/c 2 ) 6 &° oc ( GM/c 2 ) 2 R 4 5 6 . Similarly, matching the c~ 7 scaling 
with the x~ 5f B q terms in Table IIIII returns an expected scaling of ( GM/c 2 ) 5f B° oc ( GM/c 2 ) 2 R 3 for < B q . Either one 
(but not both) of these estimates could be off if the suppression of post-Newtonian order observed in the case of £ q 
also occurred here. Thus, for example, if the terms appeared only at order c~ 9 in the post-Newtonian metric, 
then Q3 q would scale instead as (GM/c 2 ) 5 ©” oc ( GM/c 2 ) 3 R 2 . 

Based on this combination of definite results and educated guesswork, we obtain the expected scalings of the 
rotational-tidal Love numbers with the body’s radius R. We express this as 
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(5.38) 


where c q , f°, b q , and t° are scalefree versions of the rotational-tidal Love numbers. These are expected to be approx¬ 
imately independent of GM/(c 2 R) when the body is weakly self-gravitating, but to acquire a dependence upon this 
quantity when the internal gravity becomes strong. We recall the caveat made in the previous paragraph: our con¬ 
siderations cannot guarantee that the expected scalings of Q3 q and A° (and even 2: q ) are not altered by a suppression 
of post-Newtonian order. 
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